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Code division multiple access (CDMA) is a spread-spectrum technology, in which a code is used 
to spread the spectral content of information to achieve larger network capacity. It has been widely 
used in classical communication and wireless networks. Here, we extend CDMA to the quantum 
case, where many pairs of nodes can transmit quantum information through a shared channel. This 
can be achieved by chaotic encoding of quantum information to spread its spectral content and by 
chaos synchronization, to separate different sender-receiver pairs. The proposed method is robust 
to low-frequency noise, and achieves faithful transmission with fidelities of up to 0.99. 

PACS numbers: 89.70.-a, 42.79.Sz, 62.25.Jk 



Introduction. — Quantum networks (for applications 
such as long distance communication and distributed 
computing) require nodes capable of storing and process- 
ing quantum information and connected to each other via 
photonic channels Recent achievements in quantum 
computing [3] have brought quantum networking much 
closer to realization. Quantum networks provide advan- 
tages for the transmission of quantum and classical infor- 
mation, after proper encoding into quantum states 0, 0] ■ 
However, the efficient transfer of quantum information 
among many nodes has remained problems yet to be 
studied [H, [a, @ • These become crucial for the limited- 
resource scenarios in large-scale networks. A remedy, 
inspired from classical communication theory, is to al- 
low simultaneous transmission of multiple quantum data 
streams to share the same channel or physical medium. 

In classical communication networks, multiple-access 
schemes have been effectively used to allow multi-user 
access. These multiple-access methods include: ei- 
ther time-division multiple- access (TDM A), frequency- 
division multiple-access (FDMA), or CDMA, in which 
each user is assigned a particular time slot, frequency, or 
code, respectively. In TDMA, the users share the same 
frequency at different time slots, but timing synchroniza- 
tion and delays become serious problems in large-scale 
networks. FDMA is bandwidth-limited, which limits the 
number of users in the network. Unlike TDMA and 
FDMA, CDMA utilizes the entire spectrum to encode 
the information from all the users by distinguishing dif- 
ferent users with their own unique codes. It is known as 
spread-spectrum technology. CDMA can accommodate 
more users per MHz of bandwidth, providing a capac- 
ity improvement on the order of 10"* to 10^, compared to 
TDMA and FDMA Q. Due to these advantages, CDMA 



achieves great success in commercial applications in clas- 
sical communication (e.g., more than 50% of the mobile 
phones currently use the CDMA technique). 

These multiple-access techniques should be applicable 
to quantum networks using either classical or quantum 
multiplexers. FDMA has already been used in quantum 
key distribution (QKD) networks [3] using conventional 
classical wavelength division multiplexers. An alterna- 
tive approach would be to introduce quantum multiplex- 
ers which transform the input quantum states from dif- 
ferent users into a quantum superposition state and then 
send it through a single channel [see Fig. fflja)]. Such 
quantum multiplexers have not been demonstrated yet. 
Moreover, despite its success in classical networks, to the 
best of our knowledge, CDMA has not been used in quan- 
tum networks and internet The use of CDMA in 
quantum networks holds the potential to provide advan- 
tages for multi-channel quantum information transmis- 
sion, like in classical communication networks. 

In this work, we extend the concept of CDMA to its 
quantum version. This extension is not straightforward, 
because in the quantum case, the information transmit- 
ted in the channel is formed by the superposition of quan- 
tum states sent by different nodes. This makes it difficult 
to coherently and faithfully extract information for each 
receiver node. Moreover, the quantum state transmitted 
through the channel is prone to noises, whose deteriorat- 
ing effects should be minimized. 

We show that obstacles towards quantum CDMA (Q- 
CDMA) can be overcome via chaotic encoding and chaos 
synchronization among senders and receivers, which also 
suppresses low-frequency noises. To achieve this, it is 
necessary to have a quantum channel to transmit quan- 
tum superposition state and N classical channels to 
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achieve chaotic synchronization to decode the quantum 
signals at the receivers. 
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FIG. 1: (color online) (a) Quantum information transmis- 
sion between two pairs of nodes via a single quantum chan- 
nel. Quantum states from two senders are combined to form 
a superposition state and input to the channel. At the re- 
ceiver side, they are coherently split into two and sent to the 
targeted receivers, (b) Schematic diagram of the quantum 
code-division multiple-access (Q-CDMA) network by chaotic 
synchronization. Wave packets from the sender nodes are first 
spectrally broadened by using the chaotic phase shifters CPSi 
and CPS2, and then mixed at a beamsplitter (BSi) and input 
to the channel. Individual signals are recovered at the receiver 
end with the help of CPS3 and CPS4, which are synchronized 
with CPSi and CPS2, respectively. 

Quantum CDMA network. — To present the underly- 
ing principle of our scheme, we consider the simplest case, 
where two pairs of sender and receiver nodes communi- 
cate quantum information via a single quantum channel 
[see Fig. [Ija)]. The quantum information from differ- 
ent users is encoded into quantized electromagnetic fields 
having the same frequencies. 

The schematic diagram of our strategy is shown in 
Fig. mb). The quantum information sent by the nodes 
1 and 2 is first encoded by two chaotic phase shifters 
CPSi and CPS2, whose operation can be modelled by the 
effective Hamiltonian 5i{t) a\ai^ with 8i{t) being time- 
dependent classical chaotic signals and i = 1,2. Here- 
after, we set h — 1. The two beams are combined by 
a 50 : 50 beamsplitter BSi and then transmitted via a 
quantum channel. At the end of the quantum channel, 
the quantum signal is amplified by a phase-insensitive 
linear amplifier (LA) and then divided into two branches 
by a second 50 : 50 beamsplitter BS2 and sent to nodes 3 
and 4 through two chaotic phase shifters CPS3 and CPS4. 
CPS3 and CPS4 are introduced to decode the informa- 
tion by applying the effective Hamiltonian —5j{t)a'jaj, 
with j ~ 3, 4. Amplifier gain is set to G = 4, to compen- 
sate the losses induced by the beamsplitters. Quantum 
interference of signals from different chaotic sources has 
recently been demonstrated [l^, suggesting the feasibil- 
ity of the scheme proposed here. 

The actions of the chaotic devices CPSi=i,2,3,4 induce 
the phase shifts exp [—i6i{t)], where 6i{t) — J* Si {T)dT. 
Thus, to achieve faithful transmission between the 
senders and the receivers, the effects of 5i{t) and 52{t) 



on the quantum signals should be minimized in the fields 
received by the nodes 3 and 4. Intuitively, this could 
be done by simply adjusting the system parameters such 
that 5i{t) = 5z{t) and 52{t) = 6i{t). However, such an 
approach is impractical, because any small deviation in 
the system parameters can be greatly amplified by the 
chaotic motion, making it impossible to keep two chaotic 
circuits with the same exact parameters and initial con- 
ditions. Instead, auxiliary classical channels between 
senders and the intended receivers can be used to syn- 
chronize the chaotic circuit as shown in Fig. [Ijb) . This 
classical chaotic synchronization enables to reduce the 
parameter differences between the chaotic phase shifters 
and to extract the quantum information faithfully. 

Hereafter, for the sake of simplicity, we assume that 
CPSi (CPS2) and CPS3 (CPS4) have been synchronized 
before the start of the transmission of quantum informa- 
tion, i.e., 6*1 (t) = 6*3 (t) [e2{t) = 6*4 (t)]. The whole in- 
formation transmission process can be described by the 
input-output relationship of the quantum network (see 
Eq. (8) in Ref. tl2j) 



04 = 02 + oie 



where aj^^ and aes are the creation and annihilation op- 
erators of the auxiliary vacuum fields entering the linear 
amplifier LA and the second beamsplitter BS2. For the 
pseudo- noise chaotic phase-shift 6i [t) , we should take an 
average over this broadband random signal 11 1| . which 
leads to exp (±16*^ (t)) sa y/Ml (see Eq. (3) in [1? 
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In Eq. ([2]), Ss^{uj) is the power spectrum density of the 
signal 5i{t)^ and lou and Uui are the lower and upper 
bounds of the frequency band oi5i{t), respectively. Equa- 
tion (P) can then be reduced to 
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For a chaotic signal with broadband frequency spectrum, 
the factor Mi is extremely small, and can be neglected in 
Eq. ([3]). This leads to 03 « ai and 04 « 02, implying ef- 
ficient and faithful transmission of quantum information 
from nodes 1 and 2 to nodes 3 and 4, respectively. 

The information-bearing fields ai and 02, having the 
same frequency ujc, are modulated by two different 
pseudo-noise signals, which not only broaden them in 
the frequency domain but also change the shape of their 
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wavepackets [see Fig. H^b)]. Thus, the energies of the 
fields ai and 02 are distributed within a very broad fre- 
quency span, in which the contribution of Wc is extremely 
small and impossible to extract without coherent sharp- 
ening of the ujc components. This, on the other hand, 
is possible only via chaotic synchronization which effec- 
tively eliminates the effect of the pseudo-noises on the 
fields and enables the recovery of ai {0,2) at the output 
03 (04) with almost no disturbance from 02 (ai). This is 
similar to the CDMA in classical communication. Thus, 
we name our scheme as Q-CDMA network. 

Quantum state transmission. — Let us further study 
the transmission of qubit states over the proposed Q- 
CDMA network using a concrete model. The qubit states 
{\(t>i) = VPiIsi) + Vl -Pi|ei), and |(/)2> = VP2I52) + 



y/1 — P2|e2), with pi, P2 € [0,1]}, to be transmitted 
are encoded in the dark states of two A-type three-level 
atoms; i.e., atom 1 in cavity 1 and atom 2, in cavity 2, 
as shown in Fig. [51 The qubit states are transferred to 
the cavities by Raman transitions and are transmitted 
over the Q-CDMA network. At the receiver nodes, the 
quantum states are transferred and stored in two A-type 
atoms; i.e., atom 3 in cavity 3, and atom 4 in cavity 4. 
We assume that the four coupled atom-cavity systems 
have the same parameters. Let \gi), jci), and jr^) be the 
three energy levels of atom i. As shown in Fig. [21 the 
Iffi) {fi) and je.;) — )■ jr^) transitions are coupled with 
a classical control field and a quantized cavity field with 
coupling strengths ili (t) and g. By adiabatically elimi- 
nating the highest energy level jr^), the Hamiltonian of 
the atom-cavity system can be expressed as 



(t) (c,|e,;)(5,| +c|l.9»>(e» 



(4) 



where Ci is the annihilation operator of the i-th cavity 
mode; gi (t) = gfli (t) /A is the coupling strength which 
can be tuned by the classical control field $7^ (t) ; and A is 
the atom-cavity detuning. The cavity fields Ci are related 
to the travelling fields Oi by 

ai = \/k Ci + ai_in , 02 = ^/k C2 + 02, in , 

a3,out = \/k C3 + 03, a4,out = C4 -I- 04, (5) 

where k is the decay rate of the cavity field; and ai^in, 
02, in (both in vacuum states) and 03, out, 04, out are the 
input and output fields of the whole system, respectively. 

The chaotic phase shifters CPSi=i,2,3,4 are reahzed by 
coupling the optical fields to four driven Duffing oscilla- 
tors, with damping rates 7, described by the Hamiltonian 
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fJ-^t -fit) Xt, 



(6) 



where Xi and pi are the normalized position and mo- 
mentum of the nonlinear Duffing oscillators, uj^jl-K is 
the frequency of the fundamental mode, is a nonlin- 
ear constant, and /(t) = cos (wrft) is the driving force. 
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FIG. 2; (color online) Schematic diagram of quantum state 
transmission over quantum CDMA network. 



The interaction between the field and the i-th Duffing 
oscillator is given by the Hamiltonian Hi — gt-oXialai, 
where gf_o is the coupling strength between the field and 
the oscillator. In optomechanical systems, such an inter- 
action Hamiltonian can be realized by coupling the op- 
tical field via the radiation pressure to a moving mirror 
connected to a nonlinear spring (see Fig. [2]) . Under the 
semiclassical approximation for the degrees of freedom of 
the oscillator, the interaction Hamiltonian Hi leads to a 

' N ' ^j^g field ai. To 



phase factor exp 



; /g 5f_oXi(r) dr 
simplify the discussion, we assume that all of the four 
Duffing oscillators have the same wq, /i, fd, and uJd, but 
different initial states. Finally, the chaotic synchroniza- 
tion between CPSi (CPS2) and CPS3 (CPS4) is achieved 
by coupling two Duffing oscillators by a harmonic poten- 
tial V (xijXs) — kj {xi — xs)^ /2, which can be realized, 
e.g., in mechanical systems, by coupling the two oscilla- 
tors via a linear spring with stiffness kj. Such kind of 
synchronization of mechanical oscillators have been real- 
ized in experiments (see, e.g., Ref. [l^). 

To show the efficiency of state transmission in Q- 
CDMA, let us calculate the fidelities Fi = {4'i\P3\4>i} 
and F2 = (</'2|p3|02), where ps and p4 are the quan- 
tum states received by atom 3 and atom 4, and j^i) = 
Vpo\9i) + \/l -Pojei) and |(/)2) = y/1 -P0I52) + VPo|e2) 
are the two quantum states to be transmitted. By de- 
signing the control parameters gi (t), using the method 
given in Ref. [3| , we find that the fidelities Fi and F2 are 
independent of the transmitted quantum states param- 
eterized by pq and just depend on the correction factor 
M — Ml — M2 ^ 1 which can be approximated as 
Fi — F2 1 — M, by neglecting the terms (see 
Eq. (33) in When the Duffing oscillator enters 

the chaotic regime, we have M f« 0, leading to fideli- 
ties Fi, F2 ^ 1, which means that the qubit states can 
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be faithfully transmitted over the Q-CDMA network. 

Numerical simulations. — Here we show the feasibil- 
ity of the Q-CDMA scheme using numerical simulations 
with the system parameters lo^i/luq = 5, .gf-o/^o = 0.03, 
fi/iOQ = 0.25, j/uJo = 0.05, ki/uiQ = 0.1, and po — 0.6. In 
Fig. EJa), it is seen that there are three distinct regions 
representing how the chaotic motion affects the fidelity of 
the quantum state transmission. In the periodic regime 
(0 < /d/wo < 15), both Fi and F2 experience slight in- 
creases with increasing fd/(^o, with 0.4 < Fi < 0.5 and 

0. 6 < F2 < 0.64. At /d/t^o = 15, the soft-chaotic regime 
of the Duffing oscillator start, indicated by a positive 
Lyapunov exponential A ~ 0.038 and a sudden jump in 
fidelities. In this regime, delineated by 15 < fd/oJo < 33, 
both Fi and F2 are still below 0.7. The dynamics of 
the Duffing oscillator enters the hard-chaos regime at 
fd/ujQ ~ 33, where both Fi and F2 jumps suddenly to 

1, which corresponds to an almost 100% faithful state 
transmission. 
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FIG. 3: (color online) (a) Fidelities Fi and F2 versus the 
strength fd of the driving force acting on the DufHng oscillator 
with Po = 0.6, and r — 2-k jioo as the unit of time, (b) 
and their average (F\ -\- F2) /2 versus po in the hard-chaotic 
region, with fd/ijJo = 36. The average fidelity is maximized 
at Po = 0.5, which corresponds to = \4'2)- 

In Fig. ISlJb), we plot the trajectories of Fi and F2 as 
a function of po in the hard-chaotic regime fdl^a — 36, 
corresponding to M « 0.0103. It is seen that Fi and F2 
are very close to 1 — 7\f « 0.9897 and almost constant 
regardless of the value of pq. There are small deviations 
from 1 — M, because here terms are not neglected. 
The average fidelity F = [Fi -I- F2) /2 is maximum at 
Po = 1/2, which corresponds to an equally- weighted su- 
perposition of the quantum states \4>i) and \4>2)- In such 
a case, the crosstalk between the channels becomes min- 
imum, inducing only a very slight disturbance on these 
indistinguishable states. 

Robustness to low-frequency noise. — Our scheme is ro- 
bust to low-frequency noises in the quantum channel be- 
tween BSi and the LA of Fig. [5] To show this, we express 
the quantum field transmitted in the channel by the an- 
nihilation operator c — c + n, where c — e*^^ai/\/2 



€'"^02/^/2 is the annihilation operator of the informa- 



i02, 



tion bearing field, n = Ec^g^o.^^.o] (^"^ + h.c.) is 
the noise term, denotes the amplitude of the noise 
component with frequency w, and WiOi '-^uo are the lower 



and upper bounds of the frequency band of the noise n, 
respectively. Then the fields entering the cavities 3 and 4 
of Fig. [5] are expressed as = + n exp [i9i (t)] , with 
given in Eq. ^ and i = 3,4. Since the noise n is mainly 
distributed in the low-frequency regime, the phase shift 
6i (t) can be divided into a low-frequency part On (t) and 
a high-frequency part O^i (t), where O^i {t) is the com- 
bination of those components much faster than those of 
n (t), i.e., the components with frequencies uj G [uiu,ujui] 
and iuii ^ u!uo. By averaging over the high-frequency 
components, we have Oi — Oi + \/ Mine^'**-*\ where 



Mi = exp 



-TT J^"' duSfi^iuj) jiJ^ . Since Mi is still very 

small, we have Oi ~ a,;, which means that the noise ef- 
fects are significantly reduced by chaotic synchronization. 
Additional calculations show that the noise n leads to 
high-order correction terms like M\M and M2M which 
can be neglected, without affecting F\ and F2. Thus, we 
conclude that our scheme is robust against low-frequency 
noises in the quantum channel. 

Conclusions. — We have introduced a Q-CDMA net- 
work based on chaotic synchronization where quantum 
information can be faithfully transmitted with fidelities 
as high as 0.99 between multiple pairs of nodes shar- 
ing a single quantum channel. The proposed quantum 
multiple-access network is robust against low-frequency 
noises. Chaotic phase shifters and their synchroniza- 
tion required for the proposed Q-CDMA can be imple- 
mented in various systems, such as optomechanical and 
solid-state quantum systems. In particular, whispering- 
gallery- mode optical resonators (e.g., microtoroid, mi- 
crodisk, microring, etc) are ideal platforms. Chaotic be- 
havior in an optical- fiber-coupled microtoroid resonator 
has been reported in Ref. [ij], and the synchronization 
of mechanical modes in two spatially-separated micror- 
ing resonators coupled via waveguides has been proposed 
in Ref. [lH, suggesting that our proposal is not far from 
practical realization. Owing to the rapid pace of experi- 
mental and theoretical developments, we believe that our 
proposal will soon be able to pave the way for long dis- 
tance Q-CDMA networks, and will give new perspectives 
for the optimization of quantum networks. 
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This supplementary material presents derivations on 
the following topics: (1) how to average out the chaotic 
phase shift; (2) derivations of the input-output relation- 
ship of the quantum CDMA network; (3) derivations 
about the quantum state transmission over the quantum 
CDMA network; and (4) analysis of the chaotic synchro- 
nization of the Duffing oscillators. 



Since the chaotic signal 6i(t) is mainly distributed in the 
high-frequency regime, we have Aia uJia- Thus, from 
the approximation Jo (a:) ~ 1 — /A, log (1 -I- a;) « x for 
X <^ 1, it can be calculated that 
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AVERAGING OVER THE CHAOTIC PHASE 
SHIFT 



Since the chaotic signal 5i{t) can be expressed as a 
combination of many high-frequency components, i.e.. 
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^iit) = ^ia cos (uJiat -|- . 
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where 



where Aia, 



bia are the amplitude, frequency, and 
phase of each component, we can obtain 

9i(t) = / di{T) dr sin {ujigt -\- (jj^g) ■ 
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From Eqs. ([2]) and ([3|), we can obtained the equation 



Using the Fourier-Bessel series identity 
exp(ia;siny) = J„(x) exp (my), with Jnix) as 

the n-th Bessel function of the first kind, it can be 
calculated 



exp [-i9^ (i)] = VM- 
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Let us take average over the "random" phase 9i it) [3| , 
the components related to the frequencies uiia should 
be looked as fast-oscillating terms and thus can be av- 
eraged out. The above treatment corresponds to av- 
eraging out the components that are far off resonance 
with the information-bearing field, and keep the near- 
resonance components. Only the lowest-frequency terms, 
with Ug = 0, dominate the dynamical evolution. Thus, 
we have 

exp [~i9i{t)] =W h ^ 



INPUT-OUTPUT RELATIONSHIP OF THE 
QUANTUM CDMA NETWORK 

Let us calculate the input-output relationship of the 
quantum CDMA network shown in Fig. [T] 

As shown in Fig. [1] the input-output relationships of 
the chaotic phase shifters CPSi=i^2,3,4 can be expressed 
as 
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and the input-output relationships of the two beam split- 
ters BSi and BS2 and the linear quantum amplifier "LA" 
can be expressed as 
BSi: 
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FIG. 1: (color online) Schematic diagram of the quantum 
CDMA network considered here. The black dashed lines de- 
note the desired chaotic synchronization channel. The red 
lines show the quantum optical channels. "LA" refers to lin- 
ear amplifier. "BS" refers to beamsplitter. "CPS" denotes 
chaotic phase shifter. 



with frequency LOd and coupling strength g. By adiabati- 
cally eliminating the energy level \ri) (see, e.g., Ref. [5|), 
we can obtain the following simplified Hamiltonian of the 
i-th coupled atom-cavity system 

H^ = ft {t) (c^\e,){g,\ + c||ft)(e,|) , (10) 

where is the annihilation operator of the i-th cav- 
ity mode; ft {t) = gVli (<) / (2A) is the reduced coupHng 
strength between |ft) and jci); A = — Wq is the de- 
tuning of the frequency between the cavity mode and the 
atom; and = — Eg = Er — E^ is the energy dif- 
ference between the ground states \ei) (|ft)) and excited 
state \ri) of the three-level atom. Here, to simplify our 
discussions, we have omitted the frequency-shift terms in- 
duced by the atom-cavity interaction, which can be easily 
compensated by tuning the frequency ujd and introducing 
an additional controllable phase-shift term by the clas- 
sical control field acting on the z-th atom (see Ref. [^). 
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= VSal + 2aLA, (7) 
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03 = -^07 -f -^flBS, 04 = -^07 - -^flBS- (8) 

From Eqs. ([5][5]), we can obtain the total input-output 
relationship of the quantum network 



a, = a2 + aie'('^-«^) + ^e*«^a[A-^e^'^aBs.(9) 



Here 61 and 62 come from unsynchronized chaotic signals 
and thus can be looked as independent "noises" . 



QUANTUM STATE TRANSMISSION OVER THE 
QUANTUM CDMA NETWORK 

Let us now derive the dynamical equation of the quan- 
tum CDMA network used to transmit quantum states. 
Here, we use the Cirac-Zoller-Kimble-Mabuchi strategy 
given in Ref. 0] to transmit quantum states. The quan- 
tum states l^i) = V^lfi) + VI - ^o|ei) and |02) 
y/1 — ro\g2) + v^|e2) are encoded in the dark states of 
two A-shaped three- level atoms in cavity 1 and cavity 2. 
We want to send them to the other two A-shaped three- 
level atoms in the cavity 3 and cavity 4 (see Fig. [2]) . Let 
Ift), jci), \ri) be the three energy levels of the i-th atom. 
The Ift) — > \ri) transition of the atom i is coupled with 
a classical control field with coupling strength fi^ (t) and 
the jci) — \ri) transition is coupled with the cavity mode 

















1 1 n-^ 


Synchronization 


ail 




Cavity 2 

FIG. 2: (color online) Schematic diagram of quantum state 
transmission over a quantum CDMA network. The broom- 
shaped or shovel-shaped purple symbols denote photon de- 
tectors. The red arrow inside each (green) cavity denotes the 
classical driving field with amplitude Qi (t) (i = 1,2,3,4). 
The green circles denote A-type three-level atoms. 

The complete Hamiltonian of the whole system shown 
in Fig. [3] is 



1 

(^ciaian - c\aiM^ + i^/k [c2a\^^ - c\a2. 



i=l 



(11) 



where k is the decay rate of the cavity mode introduced 
by the external input fields ai^in, 02, in, 0,3, and 04. From 
the input-output theory, the output fields ai and a2 of 
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the cavity 1 and cavity 2 can be expressed as 

ai = y/KCi + aia„, 02 = \/kc2 + a2^n- (12) 

Recall that the input-output relationship of the quantum 
CDMA network can be written as 
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Substituting Eqs. (HH) and ((131) into Eq. ([TTl), the total 
Hamiltonian can be re-expressed as 
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We assume that aun, a2,in, ^la, and obs are all in vac- 
uum states. The present model is a particular example of 
the cascaded quantum system which can be described by 
the input-output formalism of such a system (see, e.g., 
Ref. [1, Hi)- In the language of quantum trajectories Q, 
the dynamics of the total system can be described by a 
non-Hermitian effective Hamiltonian 

4 4 
Hcs = '^Hi-iK ^ clcj + 4 (ci + MC2) 

1=1 L i=l 

+4(C2+Afci)] (15) 
and two quantum-jump terms with Lindblad operators 

il = \/k(ci -I- C3) , L2 ^ {C2 + C4) . (16) 

Here, we have omitted the dissipation terms proportional 
to the correction factor M <^ I. Just like in Ref. 
the following dark state condition should be satisfied to 
achieve coherent-state transmission 

(01+03)1^-) = (C2+C4)|V') =0, (17) 

where is the state of the total system. 



If the correction factor M 
nian can be expressed as 



0, the effective Hamilto- 
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The quantum channel composed of node 1, node 3, and 
the corresponding cavity modes represented by ci, C3 is 
separated from the quantum channel composed of node 
2, node 4, and the cavity modes C2, C4. The state of the 
total system can then be assumed to be in the separable 
state 



\ij {t}) = 1^13 (t)) <E> \^24 it)), 
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where 

IV'13) 

IV'24) 



{vPil5i53)|0i03) + v^l -Pi [ai|ei53)|0i03) 
+a3 15163) IO1O3) -f /3i 13133) I I1O3) 

+/33|5l53)|0ll3)]}, 

|v^l5234)|0204) + VT~p^[a2 16254)10204) 
+0415264)10204) + /32 15254)11204) 

+/34|5254)|02l4)]}. 

We can find the following evolution equations by com- 
paring the coefficients of the Schrodinger equation 

ai = -9i{t) f3i, ^3 = -53 (0^3, 
Pi ^ gi{t) ai + k(33, ^3 = g3{t)a3- k(3i, 

"2 = -51 (i)/32, d4 = -53(0^4, 
/32 = 51 (i) "2 + '«/34, Pi= 9'i{t)cfi- K.f32- (20) 

Here, we have set 32 {t) = 5i if) and 34 {t) = 33 {t), which 
is reasonable because we have assumed that the four cou- 
pled atom-cavity systems are with the same system pa- 
rameters. To solve Eq. ((20)) . we can design the controls 
gi (t) {i = 1,2,3,4) by the Cirac-Zoller-Kimble-Mabuchi 
method (see Ref. [i]) 

53 it) = 91 i-t) , Vi 
g^(^^t)^-^^-PM^)+9At)aiit)^ 
as [t) 

where ai, /3i3 can be obtained by solving the equations 

ai - 51 (i) /5i3/\/2, 

/3i3 = -Kl3i3-V2gi{t)ai, (22) 
with initial conditions 
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/3i3(0) = 
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and (t) can be obtained by 



«3 it) = ^Jl-aj (0-/3?3 it). 



(24) 



By such a design of the control pulses, we can realize a 
perfect state transmission from node 1 to node 3, and 
node 2 to node 4, such that 



ai (— oo) = as (+oo) = 1, 
a2 oo) = a4 (+oo) — 1. 



(25) 



For non-zero correction factor Af <C 1, the effective 
Hamiltonian can be expressed as 



inM (^cl 
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Since M ^ 1, we can assume that the system state still 
approximately resides in a separable state, but should in- 
clude some other two-photon terms induced by the cross 



term M I 



3C2 



Ici) in iJoff (see Fig. EJ. 




FIG. 3: (color online) Intersection of the two quantum chan- 
nels when the correction factor M 7^ 0. 

The system state can then be expressed as 



I Vpi 1.91.93) IO1O3) + -Pi [ai|eig3)|0i03) 
+a3|5ie3)|0i03) -H/3i|<7i53)|li03) 

+/33|5lff3)|0ll3)]} {^/pi|<?254)|0204) 

+ y^l -p2 [a2 162.94)10204) + a4|g2e4)|0204) 
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(27) 



By comparing the coefficients of the Schrodinger equation 
and letting gi (t) = 32 (t) and 53 (t) = (t), we can find 



the following dynamical equations 

"1 = -51 03 = -53 (i)/33, 

A = 51 (t) "1 + K/33, = 53 (i) "3 - K/3i, 
"2 = -51 (i) /32, 54 = -53 (0/34, 

/32 = 51 (t) "2 + K^4, ^4 = .93 (t) "4 - k/32. 



(28) 
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"3, 
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v/p2(l-pi)/VPi(l-P2)l /33- (29) 



/34 = ^4-l-Af 



Comparing Eqs. (pHl) and (pS)) . we can see that the dy- 
namical equations for a^, /3i, i = 1,2,3,4 and the equa- 
tions for ai, a2, as, a4, /3i, /32, /Ss, /34 are the same. We 
can design the controls gi (t) [i = 1, 2, 3, 4) by simply re- 
placing as, a4, /3s, Pi, by as, 04, /3s, and Pi in Eqs. ([2T|) . 
(I22|), and dSSl). Then we have 



(30) 

0(Af2), 
0(Af2), (31) 



as (+00) = a4 (+00) = I, 
which means that 



as (+00) = 1 - Af 
a4 (+00) = 1 - A/ 



yjpi (1 -P2) 



\Jvi (1 -Vi) 
\/vi (1 -Vi) 



\Jv\ (1 -V2) 



where O(Af^) represents higher-order terms of M under 
the condition that A/ <C 1 . It can be further checked that 
^5 (-1-00) = /36 (+00) = 0. To solve as (-I-00), ag (-I-00), 
we can see that the reduced state of the system composed 
of atom 3 and atom 4 can be written as 

P34 (0 (\/Pi|.93> + -Pias (+00) les)^ 



( 



'P2I54) + -P2a4 (+00) |e4) 
(\/Pi(.93| + v^l -Pias (+00) (esl) 
^\/P2(54| + a/1 -P2a4 (+00) (64]^ 
(+00) 16354) (es54| 

-ag (+00) |5364)(5364|- 



(32) 



From tr[ps4 (i)] = 1, we have 

ag (+00) + ag (+00) 
- 1- [pi + (l- Pi)a2 (+(»)] [p2 + (l- 

= 2A/ VpiP2 (1 - Pi) (1 - P2) ( — + — 

VPi P2 



P2)a4 (+00)] 
+ 0(Af2). 
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Further, from the symmetry of the system Hamiltonian 
(I26p . we have (+00) = ctg (+00), which means that 

al (+00) = Q!g (+00) 

= M VpiP2 (1 - Pi) (1 - P2) (— + —)+ 0(M2). 

\Pi P2J 



From Eq. (1521) , it can be shown that the reduced states of 
atom 3 and atom 4 in the long-time hmit can be expressed 
respectively as 

P3(+Oo) = [p2 + (1 - P2) "I (+00)] 



Pilgs) + y/l-pias (+00) lea) 
Pi (53 1 + Vl -pias (+00) (63 1 



/94 (, + CX)J 



(+00) |e3)(e3| + al (+00) \g3}{g3\, 
[pi + {l-pi)al (+^)] 

P2I54) + -P2Q;4 (+00) |e4) 



^2(54! + Vl -P2a4 (+00) (64 1 

+al (+00) |e4)(e4| + (+00) |54)(g4|- 

Thus, we can calculate the final fidelities of the state 
transmission through the two channels 
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(VPilffa) + Vl -Pi|e3)) 
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[pi + (1 -_Pi)q!3 (+00)]^ 
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\Pl P'. 

I- M^pip2{l-Pl){l-P2) ( — + - 

VPi P'. 

+0{M^) 

(VP2(54| + ^1-^2(64!) P4 (+00) 

(^^154) + -P2|e4)) 

[P2 + (1 - J32) "4 (+00)]^ 
+M ^JpiP2 (1 -J3l) (1 -P2) ( — + — 

1 



1 -Mv/Pip2(l-Pl)(l-P2) ( — + — 

VPl P2 

+0{M^). 

If the transmitted quantum states are 
|0i> = Vpo\9i) + -po|ei) 

and 



(33) 



i.e., pi = po, P2 = 1 ^ Po, it can be calculated from 
Eq. dMl) that 



Fi = F2 = 1 - M + 0{M'^). 



(34) 



CHAOTIC SYNCHRONIZATION OF DUFFING 
OSCILLATORS 

We will analyze the chaotic synchronization of the 
two driven Duffing oscillators shown in Fig. |4] induced 
by the harmonic interaction potential V{xi^X3) = 
kj {xi — x^)^ /2 , and give an estimate of the bound of 
the interaction strength kj. 
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FIG. 4: (color online) Schematic diagram of the chaotic syn- 
chronization realized by the moving mirrors. 

The dynamics of the two coupled Duffing oscillators is 
described by 

7 

xi = wopi - -xi, 

pi = -uJoXi+4fiXi + f{t)-kj{xi-X3)- -pi, 

1 

X3 = ^oPs - 2^3, 

3 7 

P3 = -woa;3 + 4^X3 + /(t) - fc/ (0:3 - xi) - -p3(35) 

From these, we can obtain the equations of the differences 
of the positions and momenta of the two oscillators, = 

Xi — X3 and Bp — Pi ~ P3, as 

■ - _ 2 



where 



Bp ^ (-cjo + 4:^X - 2ki) ex - -e^ 



X = Xi + X1X3 + x^. 



h) = Vl -po\92) + VPo\e2), 



Since the eigen-equation of the linear matrix of the above 
equation is 
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it can be verified that the chaotic synchronization, i.e., 
63;, Bp — 0, occurs when 
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ki > AfiB - Wo, 

4wo 



(36) 



where B is the upper bound of \X\, i.e., \X\ < B. Equa- 
tion ((36|) shows that the two chaotic oscillators synchro- 
nize only when the strength fc/ of the harmonic interac- 
tion potential V (xi, xa) is strong enough. 

To estimate the lower bound of fc/ to achieve chaotic 
synchronization, let us give an estimation of B. Since 

1^1 < \xi\'' + \xi\\xz\ + \xz\'', 

we just need to estimate the upper bounds of \xi\ and 
1 2:3 1 . Assume that the two Duffing oscillators have syn- 
chronized, then we have (xi — 0:3) — >■ 0. Thus we only 
need to estimate the bound of xi, and omit the cross 
term — fc/ (xi — X3) in Eq. (j35p to obtain the equations 
of Xi and pi as 



Xl 
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Pi = — wqXi + 4^Xj 



(37) 

To estimate the bound of xi, we introduce the har- 
monic balance approach, by which we can omit the high- 
frequency components and assume that 



Xl — A cos {uJdt - 



(38) 



Substituting Eq. ((38)) into Eq. (|37| . we can obtain the 
following equations by comparing the coefficients of the 
terms cos {udt -\- 4>) and sin {uJdt + </*) 
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3^A3 



fd cos ( 



Eliminating the unknown parameter 0, we can obtain 
the algebraic equation to solve the bound A of |xi| 



A' 



1 
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UJQ 



n 2 



LdQ + 3p,A^ 



(39) 



Then, we can verify with simple calculations that the 
two Duffing oscillators synchronize only if the coupling 
strength kj between them is large enough such that 



ki > UnA^ 
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4^0 



(40) 



where the parameter A is solved numerically from 
Eq. 1^. 

We now show some numerical results, using the system 
parameters 

(wd, 5f _o, ^J., 7, ki) = (5a;o, 0.03tJo, 0.25a;o, O.OSwq, O.Iwq) ■ 
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FIG. 5: (color online) Chaotic synchronization: (a) the blue 
curve shows the difference {xi—X4,) (the unsynchronized case) 
and (b) the black curve shows the difference {xi — x^), with 
ki = O.lwo (the synchronized case). 



In Fig. m we show the chaotic synchronization induced 
by the harmonic interaction potential y (xi,X3). Here, 
we set /(i/wo = 36, where the dynamics of CPSi=i_2,3.4 
(the Duffing oscillators) are all in the chaotic regime. In 
Fig. Efa), the normalized position xi of CPSi and X4 
of CPS4 are unsynchronized in the long-time limit due 
to the lack of the harmonic interaction. In this case, the 
small difference in the initial condition is amplified in the 
long-time limit due to the sensitivity to the initial con- 
dition of the chaotic motion. In Fig.jSjb), the difference 
between the position xi of CPSi and X3 of CPS3 tends to 
zero, which means that xi and X3 are synchronized under 
the harmonic interaction potential y(xi,X3). From the 
parameters chosen in our simulation, it can be estimated 
from Eq. ^ that kj > O.OSujq. Thus, the two Duffing 
oscillators synchronize with kj — OAujq in our simulation 
results shown in Fig. [5Kb). 
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